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Summary. The comparison between observations of galaxy clusters thermo-dyna- 
mical properties and theoretical predictions suggests that non-gravitational heating 
needs to be added into the models. We implement an internally self-consistent 
heating scheme into GADGET-2 for the third (and fourth) run of the Millennium 
gas project (Pearce et al. in preparation), a set of four hydrodynamical cosmological 
simulations with = 2 x (5 x 10^) particles and with the same volume (L = 
500/i~^Mpc) and structures as the the N-body Millennium Simulation (Springel et 
al. 2005). Our aim is to reproduce the observed thermo-dynamical properties of 
galaxy clusters. 

1 Model 

The large dynamical range that characterises cosmological simulations and 
the fact that the physics of heating mechanisms like AGN feedback, galactic 
winds and conduction has typical scales much smaller that those describing 
galaxy clusters (few kpc versus Mpc), make their implementation challenging. 
In addition the Millennium gas runs do not have a high enough resolution 
to properly model these phenomena {nigas — 3.12 x 10^/i~^Mo, softening 
25/i~^kpc) and therefore we need to adopt a relatively simple heating scheme. 
Our model seeks to improve on the simple preheating scheme that was im- 
plemented in the second run of the Millennium gas project. 
We implement a self-regulated star formation plus feedback procedure by 
selecting gas particles to convert into stars by imposing a temperature plus 
over-density threshold (in units of critical density) . We choose an over-density 
rather that a physical density threshold. We then inject an energy Einj to 
the neighbours of the new star, weighted in distance by the standard SPH 
smoothing kernel. 

We choose a temperature threshold of IQ^K in both models and density 
threshold, pthr, to be: 

~ Pthr = 200: mimics the energy deposited by supernovae in galaxy cluster 

outskirts, low mass objects and at early times 
~ Pthr = 2500: mimics the effect of star formation in the centre of galaxy 

clusters 
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We then tune Einj in order to reproduce the observed luminosity-temperature 
relation of galaxy clusters at z = 0. 




Fig. 1. Left panel: bolometric luminosity versus mass weighted temperature for 
groups and clusters in a test box of size L = 125ft~^Mpc. Black dots are for the 
preheating run, red diamonds for pthr = 200. Observations from Ponman et al. 
1996; Helsdon et al. 2000 and Novicki et al. 2002. Right panel: hot gas (shaded 
area) and stellar fraction (lines) for the preheating run (black lines) and pthr = 200 
(red lines). 



2 Results 

The amount of energy that we inject, Einj, is our free parameter and it will be 
different for the two runs, which vary the density of gas into which the energy 
is injected. From a test run with L = 250/i~^Mpc we found that choosing 
Einj = l4:keV for pthr = 200 we manage to reproduce the observed L — T 
relation for group size objects as well as clusters in the range of 1 — 6keV. 
Comparing this run with the preheating model (similar to the model proposed 
by Kay et al. 2004), we find that the two runs have a slightly different nor- 
malisation but still in agreement with observations (Figure 1). We also notice 
that the L — T relation has too little scatter in the preheating run, compared 
to the observations, while the current model appears far more reasonable. 
This difference is at least partially due to the difference in the temperature 
structure: in the preheating run the gas is brought to such a high temperature 
at early times that no cool gas is found at redshift zero and the haloes are 
characterised by a very smooth gas distribution. On the contrary cool cores 
can form in our current model and the clusters often show multiple-structure 
in temperature maps and offsets between the emission peak and the temper- 
ature peak (Figure 2). 
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Fig. 2. Mass weighted temperature maps for three clusters. The panel size is four 
times the virial radius. Superimposed surface brightness contours. 



In the second panel of Figure 1 we show the hot gas fraction and the 
star fraction as a function of cluster mass. Tuning the preheating run and 
pthr — 200 to match the observed L — T relation at z = 0, results in a much 
reduced baryon fraction relative to the non-cooling case, with fi,ar < 0.13 
while the universal fraction is 0.18 for this cosmology. It is interesting to 
notice that we obtain similar trends despite the two very different heating 
mechanisms; one is an external heating (preheating run) and the other one 
an internal heating. Once the Millennium volume will be available we can use 
the simulations to test how powerful fhot{z) -hot gas fraction- and fbar{z) 
-baryonic fraction- are as cosmological probes. 

The full 500/i^^Mpc millennium run will also be suitable for time evolution 
studies. The major drawback of the preheating scheme is that the gas is 
brought to such a high adiabat at z=4, when the preheating occurs, that 
no gas can condense any more and star formation is quenched. With this 
alternative model we get a more realistic star formation history and time 
evolution. 



3 Conclusions 

We present an alternative to a simple preheating model by introducing a 
self-regulated star formation plus feedback scheme. The improvements of 
this model relative to the preheating one are: an increased scatter in the 
luminosity-temperature relation, the existence of cool core clusters and a 
more realistic time evolution. On the other hand the amount of energy re- 
quired in order to match the observational data with this scheme is quite 
high, especially if we employ a high density threshold that has been tuned to 
match the observed L — T at z = Q. We are currently exploring higher and 
lower resolution simulations to constrain the systematic effects introduced by 
our model. 
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